Ground state study of simple atoms 
within a nano-scale box 
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Ground state energies for confined hydrogen (H) and helium (He) atoms, inside a penetra- 
ble/impenetrable compartment have been calculated using Diffusion Monte Carlo (DMC) method. 
Specifically, we have investigated spherical and ellipsoidal encompassing compartments of a few 
nanometer size. The potential is held fixed at a constant value on the surface of the compartment 
and beyond. The dependence of ground state energy on the geometrical characteristics of the 
compartment as well as the potential value on its surface has been thoroughly explored. In 
addition, we have investigated the cases where the nucleus location is off the geometrical centre of 
the compartment. 
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I. 



INTRODUCTION 



In the past few years, spatially confined atoms and 
molecules have gained relentless attention 1]. This comes 
mainly from the notable difference between the ground 
state properties of such confined systems in comparison 
to their bulk states. As a result, the confining systems 
may be applied in several circumstances. At present, it is 
argued that physics of high pressure materials can be in- 
vestigated within models of confining particles. Further- 
more, extensive investigations over large classes of nano- 
structure systems such as quantum dots and artificial 
atoms are essentially related to the problem of confined 
atoms and electrons By the issue confined atoms, 
we mean the atoms that experience external potentials 
which keep themselves in a region with length scales com- 
parable to the atomic size. Certainly the boundaries of 
such regions are not fully impenetrable, and particles are 
then able to escape from the region. Generally the con- 
fining surface is not spherical. The energy spectrum of H 
and He- atom in penetrable/impenetrable spherical and 
ellipsoidal boxes have been studied by several authors 
who have employed several approaches H, 0, H, H, S H 111 ■ 
In all of these models, it is assumed that confining po- 
tential only acts up to a cut-off boundary, and do not 
disturb the atomic Hamiltonian inside the region. Under 
this simplification, the problem of confined atoms, would 
be equivalent to the problem of many-body Schrodinger 
equation with a given boundary condition. Practically, 
finding the solution of many-body Schrodinger equation 
with specified boundary condition, is not an easy task 
even by numerics. As a simple but less accurate method 
to solve this equation, self consistent Hartree-Fock, or 
variational techniques are usually taken into account. 
Another methodology for solving many-body problems. 



is Quantum Monte Carlo (QMC) method. This method 
is among the most powerful methods in obtaining the 
ground state energy of a many-body system, such as 
two-electron atoms constrained in spherical impenetra- 
ble boxes p^ . In this article, we have used a familiar 
variant of QMC technique the so-called diffusion QMC 
method to obtain the ground state energy of the 

simplest confined atoms i.e.; H and He. Our motivation 
arises from the fact that in real confined atoms, the lo- 
cation of the nucleus does not exactly coincides with the 
centre of boxes. Since the analytical approaches fail in 
these cases, employing numerical techniques is inevitable. 
In order to take into account the effects of size, shape and 
the potential of the boundary on the ground state energy 
of confined H and He atoms, we have examined two ge- 
ometries namely spherical, and ellipsoidal compartments 
with different sizes and potentials at their boundaries. 
It is worthwhile to note that the problem of H-atom in- 
side a sphere or cylinder, as a model for H absorption in 
fuUcrene or carbon nanotube |T3] , has found appli- 

cations in the condensed matter physics. Moreover, the 
problem of confining He, as the simplest problem com- 
bining repulsive electron-electron interaction effect with 
the repulsion resulting from confinement, could be of in- 
terest from theoretical point of view. For example in the 
finite temperature we expect an ellipsoidal shape for Ceo 
molecule instead of a spherical geometry and in this case 
one can model this molecule by a prolat-shaped cage ■ 

Our new findings can be summarized as follows. We 
have proposed a set of trial wave functions for H- atom 
and He -atom in a penterable/impenterable spherical and 
ellipsoidal box and computed the effect of off ccntrality 
on the ground state energy for various radii and pene- 
trability. We recall that for these types of atoms, the 
number of electrons are less than three. Under this con- 
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dition, it would not be vital to use the fixed node ap- 
proximation as is used for fermionic systems [Tst . We 
first try to solve the centrally symmetric problems, and 
then will consider generalizations which involve the ef- 
fect of displacing the nucleus off the geometrical centre 
of the compartments. The paper has the following orga- 
nization: In sec. II we briefly describe our used model. 
Sec. Ill reviews the basic ingredients of DMC method 
and the specific algorithm used in the paper. In Sec. IV 
we have introduced the used trial wavefunctions and sec. 

V includes the main results. The paper is ended in sec. 

VI by some concluding remarks and discussions. 

II. THE MODEL 

Consider a simple atom subjected to a constant exter- 
nal potential Vq which acts beyond the exterior region of 
an encompassing compartment. The time-independent 
Schrodinger equation for such H-atom is: 

HtpiK) ^ Ei}(R). (1) 

where 

i?(R) = T(R) + VeniTi) + VeeCR) + 14.t(R), (2) 

here R = (r, x) and r, x denote the coordinates of the 
electrons and the nucleus respectively. T is the kinetic 
energy term for electrons. The nucleus is assumed immo- 
bile due to its large mass compared to electron's mass. 
Ven is the Coulombic interaction potential between the 
electrons and the nucleus and Vext is the external ap- 
plied potential. Vee is the electron-electron interaction 
for many electron systems which is zero beyond bound- 
aries. The external potential is zero inside the compart- 
ment and constant Vq on its surface and beyond. 



■0s (R) which is a functional of the potential T^(R) which 
is identical to the ground state wave function of the orig- 
inal quantum mechanical problem. In many situations, 
'(/'(R) shows divergent to infinity or convergent to zero 
types of behaviour. For examples if — y(R) is positive, 
as in the hydrogen atom with potential V{r) — 
Evidently ipsi^) diverges because after a sufficient time, 
the number of particle in the ensemble diverge due to 
creation process. Now to remedy the problem, one adds 
a trial constant Et into the imaginary time Schrodinger 
equation as follows: 

f^^^(^ = D\7li,{R, t) + {Et - V(R))V'(R, r). (4) 

In the limit r — > cx3, -(/'(R, t) will only tends to a non 
zero and finite value provided Et coincides with the sys- 
tem ground state Eq. To make this approach applica- 
ble, we can exploit different methods of solving diffusion- 
like equation such as random walk and Green function. 
Although a significant advantage of the diffusion Monte 
Carlo is that one does not need to specify a trial wave 
function, to speed up the implementation of the method 
and to avoid errors arising in the problems containing 
singularities in the potential, it would be helpful to use 
a guiding trial wavefunction. More specifically, in the 
positions of singularities, where T^(R) diverges, the cre- 
ation/annhilation rate ET — V{Ti) becomes tremendously 
large and can cause undesirable fluctuation in the num- 
ber of diffusive particles which makes the algorithm un- 
stable. To smooth this singularities, one uses a time 
independent trial wavefunction ■(/'^(R)- To implement 
this idea, we introduce a function /o(R, r) defined as 
p(R, r) = '0r(R)V'(R-i ''')• Substituting tACR, t) in terms 
of p in equation (4) gives the following equation: 



III. DIFFUSION MONTE CARLO METHOD 
FOR SOLVING SCHRODINGER EQUATION 

In this section we briefly outline the main ingredients 
of the Diffusion Monte Carlo scheme. The essence 
of DMC lies in mapping the Schrodinger equation in 
imaginary time onto a diffusion-like equation [ll|, [isj . 
Through an analytic continuation of time < — > — it to 
imaginary values, the Schrodinger equation will appear 
in the following form: 



' dr 



W^VCR,^) -F(R)V'(R,t) 



(3) 



In which the constant D is This equation describes 
the diffusion of an ensemble of interacting particles 
subject to death and/or growth processes depending 
on the sign of V^(R). The rate of death or growth i.e.; 
— V^(R) is not constant but rather space-dependent. In 
the long time limit, ^(R, r) tends to a unique value 



' dr 



Wr[Vr-F(R)]p(R, r) 
{El{B.)~Et)p{B.,t). 



(5) 



The term Fili) = ^^^[ff^ is called Fokker- Planck 
force and corresponds to force function which drifts 
the walkers away from regions where |?/;t(R)P is 
small. The local energy function i?i(R) is defined as 

~^^t(r/^'' ~^ ^(R)- appropriate choice of the trial 
wavefunction '0t(R) gives rise to non singular _Bl(R). 
The method we have taken into account is based on 
the Green's function approach to equation (5). We 
recall that the Green's function of equation (5) can 
not be obtained exactly. However, the problem can be 
treated using an approximate Green's function for short 
time evolution. The appropriate approximated Green's 
function is divided into two terms [i8[ which are given 
below: 
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Gz,(R,R', Ar) = ^.^^^fj/. ^ + 0(Ar)^(6) 



(47rl? At)3^/2 



for diffusive term and 



G.(R-,R, Ar) = exp[-Ar( ^^(^) +^^(^^) - Et)1 (7) 

for the branching term. This is the part of the 
Green's function responsible for the stochastic cre- 
ation/annihilation process. To improve the error of the 
order (Ar)^, we have used the Metropolis procedure 
which will be explained shortly. Throughout the paper, 
we have used the following atomic units h = m = e = 1. 
We note that the atomic unit of length is the Bohr radius 
= 0.529 A : the time unit is = 2.419 x IQ-^ fs 



and the energy unit is Hartree i.e. 



me 



= 27.25 eV. 



and with the probability 1 — {q — [q]), [q] replicas are 
replaced at the updated location of the particle to the 
ensemble of particles. Note that if q is less than one, 
with probability 1 — q the particle annihilates. After 
performing the above diffusion/branching processes for 
all the particles in the ensemble, and hence updating 
the ensemble number TV, the final MC step consists of 
updating the value of Et- Among various choices of 
updating rule, we have chosen the following one: 



Ex =< E]^ > 



N-No 
' No At 



(10) 



Where Nq refer to the initial number of walkers, N de- 
notes the updated number of walkers. The average is an 
ensemble average. This adjusting of Et is essential to 
avoid large fluctuation in the walkers number. After a 
sufficient MC steps, the time series Et reaches a steady 
state. The ground state of the system is obtained by 
averaging Et over many MC steps. 



A. DMC steps 



IV. TRIAL WAVEFUNCTIONS 



In what follows we give the algorithmic steps for 
implementation of the method described above. Firstly, 
we initialize the system. A number Nq of non interacting 
particles (random walkers) are distributed in space. 
It is preferable to distribute them according to an 
approximate guess of the ground state wavefunction. 
Then the position and number of walkers are updated 
according to the following stochastic rules. These rules 
are applied synchronously to all the walkers of the 
ensemble. During an MC step, each walker execute 
a guided discrete time random walk. To this end, we 
choose a constant time interval At between walks. The 
walk length and its direction are stochastic variables. 
Since the walk is directed, it contains a directed plus a 
random component. More precisely, the new position of 
each particle is update according to the following rule: 



R' ^ R + DArF(R) 



(8) 



In which 77 is a Gaussian random number with zero mean 
and unit variance. The term DArF(R) corresponds to 
the guided part. The attempted move is accepted if the 
Metropolis test is successful. This test introduces a ratio 
w given as below: 

^ G(R,R',Ar)p(RO _ ^(R, R', Ar)V4(R0 
^ G(R', R, Ar)p(R) G(R', R, At)V4(R) ' ^' 

The move is accepted with the probability w. The next 
stage simulates creation/annihilation processes. After 
accepting the move, we evaluate the branching factor 

q ~ 5 Let [q] denote the integer 

part of q and q— [q] its fractional part. When q is greater 
than one, with the probability q — [q], [q] + 1 replicas 



A. Spherical penetrable boundaries: Centered 
atom 

In this paper, we have used trial wavefunctions which 
are appropriately chosen associated to the spherical 
boundary condition. We note that the potential takes 
a constant value Vq on the boundary and beyond it. For 
the spherical boundary a denotes the radius of the sphere. 
The impenetrable case has already been done in |10l |. 

When the boundaries are penetrable, for H-atom, we 
have used an iterative scheme to find the appropriate 
trial wavefunction. This scheme is implemented accord- 
ing to the following steps: at first we employ the same 
unbiased method used in [ll[, for producing the radial 
wavefunction for all values of Vq and a. After obtaining 
the required data for the radial wavefunction, in the sec- 
ond step we find the best fit of the following function to 
these data. We propose the trial wave function as follows: 



(r, ^o) = Y3 expi-Yir"" - Y2r^ - Yor). 



(11) 



The fitted function ^'^(r) is in turn used as the trial 
wavefunction. We note the parameters l^'s are functions 
of Vo and a. By using the cusp condition, V't('')+V't('') = 
0, we can determine one of the above parameters such 
as Yq which is equal to one. This condition guarantees 
singularity cancelation at regions where r is small. The 
values of the fitting parameters are exhibited for some 
values of a and Vq in Table I. In this case El for r < a 
turns out to be 



ELH{r,Vo) 



{6Yir + 2Y2) - (SFir^ + 2Y2r + YoY 

1 
r 



[-SFir^ - 2Y2r -Yq + I] 
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For exterior regions where r > a the local energy should 

change to ELni'i', Vo)H hVo- Note that the local energy 

r 

is obviously finite in the region r > a. In the case of He 
atom inside a penetrable sphere the trial wave function 
is \E'^(ri)\l/^(r2) J(ri2). Here the Jastrow function has 

the form J{ri2) = e^<^^^°^T^ with ri2 as the distance 
between two electrons and a = 0.2. The cusp condition 
leads us to choose Yq = 2 . Then the local energy when 
ri < a and r2 < a is written as below 



is a complex function which in the range d < 0.5a has 
a simple linear behaviour, A = Az + B where A and B 
depend on d. Table II includes the values of A and B for 
many geometries. 

When the boundaries are penetrable and off-centrality 
is small, we have used \E'^(r',Vb) with the same func- 
tional form as in the equation (11). The corresponding 
parameters 1^'s are taken from Table I. 

In the case of small ofF-ccntrcd He atom within an im- 
penetrable spherical box we have used the following trial 
function: 



ELHe{r,VQ) = ELHiri,Vo 

2 



ELH{r2,Vo) + — 
ri2 



-E 



V,*i ViJ 



1 



J 



+ 2- 



J 



where as before J = ^(^12) and — ^'^(ri,Vb) 
according to equation (11). In the case of ri > 
a or r2 > a when one electron goes out- 
side of the sphere the local energy is written as 

ELHe{r,Vo)+ —- — 

In the above equation k refers to the index of elec- 
tron in outside region. Finally when both elec- 
trons are outside the region, i.e., ri > a and 
r2 > a the local energy should be written as 

Ti ri2 



ELHe{r,Vo) + 



B. Spherical boundeiries: Off-centred atom 

Now we discuss the off-centrality in the spherical ge- 
ometries. When the off-centrality is small and the bound- 
aries are impenetrable we have obtained, for H atom, the 
trial wavefunction by implementing a shi ft d along the 
z-direction: r' = ^/ x'^ + y"^ + {z — d^ = \J -\- [z — (£f. 



*TZ,(r') 



„(-r'-|-A(z,d)r'=) 



io(fcr). 



(12) 



where joikr) is the zero order Bcssel- function with k = 
w/a. The first term plays the role of off-centrality and the 
second term insures the vanishing of the wave function 
on the boarders. Defining / = -l-A(z,d)r .^^^ using 
cylindrical coordinates, the local energy is found to have 
the following form: 



*Ti.(r'i)*?i,(r^)J(ri2), 



(14) 



where r^' = \/xf 



Vi + (-^i — d)^ and ri2 is defined as 
the distance between two electron as before. 



C. Ellipsoidal boundaries 

In the case of ellipsoidal geometries first we recall the 

equation of a standard ellipsoid which has been pro- 
lated in the z-direction by defining the function R = 













z2 



face with the semi-major axis a and the semi-minor axis 
b. Without loosing generality we assume b> a. For small 



eccentricity i.e, e 



—5- — 1, our trial wave function for 



H atom within an impenetrable ellipsoidal boundary is 



2a' 



(15) 



These 



k should be equal to tt while /? = [2^ 

values have been determined by keeping the local energy 
finite and reaching to the same function of spherical case 
when a = b. For H atom in penetrable boundaries we 
have used (R, Vq) using the functionality analogous 
to equation (11). 

In the case of He atom in a prolated impenetrable el- 
lipsoid with small eccentricity we have Tised the follow- 
ing wave function 'I'^g(ri, i?i)^'^g(r2, i?2) J(^'i2)- For He 
atom in penetrable ellipsoidal boundaries one can replace 
the firs two terms of equation (14) by ^'y(i?i,Vb) and 
\E'^(i?i,Vb) respectively. 



ELoip, z) = 



1 

'2 
2 



^d\f/dp^ 



f 



1 rfln/ 
P 



dp 



(Pf/dz^ 



+ 



dp dz 



k cot(kr) 

r 



(13) 



where f = p + z. For keeping the local energy finite 

every where especially at r = a one needs to solve an 
ordinary differential equation for X{z,d). Such solution 



V. RESULTS 

Prior to presenting our results, it would be illustra- 
tive to discuss the numerical errors in the DMC method. 

Principally, there are two types of errors which limit 
the accuracy of most DMC calculations: (a) statisti- 
cal or sampling errors associated with the limited num- 
ber of independent sample energies used in determining 
the ground state energy and (b) the systematic errors 
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associated with finite time-step At, round-off in com- 
puting, imperfectness of random number generators etc. 
For a spacial boundary condition like the spherical boxes 
with penetrable boundaries, we know the energies value 
(steady state energy) from the analytical calculations 
([lol|). therefore, minimizing the variance of energy, can 
give us a range of time step. We have used the mini- 
mum of these time step, i.e., At=0.005. The number of 
walkers was taken iVo — 5000 in our simulation. 



A. Spherical box 

The first quantity that we have calculated is the elec- 
tronic ground state wavefunction (see Fig. 1). Due 
to spherical symmetry, the solution of the Schrodinger 
equation can be solved by separating the variables so 
that the wavefunction is represented by the product 
■0(r, 6', 0) — $(r)O(0)(^((/)). Furthermore, because the 
ground state wavefunction of centred H-atom in a spheri- 
cal geometry has radial symmetry, one can replace i/) with 
r^ir) for simplicity. For an atom confined in the centre 
of an impenetrable and very small sphere (Vq = oo and 
a=l), the wavefunction dramatically differs from that of 
free atom. Accordingly, the maximum of the wavefunc- 
tion shifts towards smaller r, which also causes larger 
slope in its curvature (see Fig. 1(a) ). On the boundary, 
the wavefunction tends smoothly to zero. Since kinetic 
energy corresponds to the Laplacian of the wavefunction, 
this enhancement of the slope consequently leads to the 
enhancement of the ground state energy. The situation 
is significantly different for large encompassing regions 
where the radius of the sphere is notably greater than one 
atomic unit. As can be inferred from Fig. 1(b), where 
a=4.0, the electronic wavefunction deviates slightly from 
the electronic wavefunction of a free H-atom. In this case 
the wavefunction tends to zero with small slope. More- 
over, the infiuence of Vq has been suppressed. 

An important feature of the electronic wavefunction 
in small compartments (when a = 1) is its dependence 
on the penetrability Vq. By decreasing the potential Vq, 
the tendency of an electron for being near the bound- 
ary increases and for pure penetrable sphere (Vb = 0.0) 
the probability of finding the electron near the boundary 
becomes maximum. In Fig. 2(a) we have plotted the 
ground state energy versus the strength of penetrability 
for several radii. Apparently, the ground state energy 
strongly depends on the penetrability when the size of 
bounded region is comparable to the volume of a free 
H-atom, i.e., a< 1. Nonetheless, for larger regions and 
larger radii, ground state energy is almost independent 
of K). 

To improve our understanding, in the set of coluored 
two dimensional figures, we have depicted the behaviour 
of electronic distribution (see Fig. 3-5) for some values of 
Vq at a = 1 in the x~ z plane. As can be seen, by increas- 
ing Vo, the electronic distribution becomes more localized 
around the nucleus and keeps its symmetry independent 



of Vq . In the case of fully penetrable boundary, electron 
has access to exterior regions. The blue regions represent 
a smaller probability for finding the electron, whereas the 
red regions represent higher probabilities. The electron 
does not penetrate to the exterior regions in deeper con- 
fining potential barrier. 



B. Off-centred atoms 

One of the most significant advantages of DMC method 
lies in its simplicity in being applied to various boundary 
conditions. This property allows us to investigate the ef- 
fect of displacement of H-atom along the z-direction on 
the ground state energy of the system. Our motivation 
in considering this geometrical case is two-fold: first, in 
real systems (such as confined atoms in a metallic bubble) 
one can not fix the atoms on the geometrical centre and 
second, due to mathematical difficulties arising from the 
analytical treatments. Figure 2(b) illustrates the results 
of this investigation for three values of Vq. The figure 
depicts the ground state energy of a H-atom confined 
in a sphere with radius a = 1 versus the nucleus dis- 
placement d (in atomic unit) from centre of the sphere. 
Generically, the off-centre effects become more consider- 
able when the displacement d exceeds nearly half of the 
radius a. Table HI contains our results for ground state 
energy. Our energy results for the case d = and impen- 
etrable boundaries are in good agreement with those in 
[To] . For the case of penetrable boundaries our results are 
in well agreement with the analytical results in For 
a deeper understanding, we exhibited the dependence of 
the energy on the displacement d in Fig. 2(c) for impen- 
etrable sphere. For larger a the deviation of the energy 
occurs in larger displacement d. 

For off centred helium atom all the results are shown 
in table IV. The general behaviour is qualitatively anal- 
ogous to the H- atom. However, due to electron-electron 
interaction, geometrical effects are suppressed compared 
to H- atom case. Our results for d = and impenetrable 
boundaries are in reasonable agreement with the results 
given in [l^] ■ It has been verified that if one uses a local 
and negative penetrability in the boundaries when the 
atom moves towards the boarder, the electron can cross 
the well depending on the well deep [lB| ■ 



C. Ellipsoidal box 

Let us now investigate the ground state energy of a H- 
and a He- atom confined within an ellipsoidal compart- 
ment. We consider a symmetric type of ellipsoid char- 
acterized by parameters a and h which are semi-major 
and semi-minor axes respectively. For simplicity, we have 
used small eccentricity, i.e, h = 1.1a and h ~ 1.5a for sev- 
eral a and Vq. We recall that nucleus is placed in the 
centre. 
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The ground state energy of H-atom confined in larger 
ellipsoidal cavity (large a and b) tends to the energy of 
the ground state of free H-atom, -0.5 , see table V. When 
we increase Vq to a higher value we observe that energy 
increases especially for smaller a. For a fixed a , larger b 
has larger energy in all cases. For an ellipsoid with small 
semi-major axis, the ground state energy will be more 
affected by the surrounding potential. 

For many-particle systems, finding an exact solution of 
the Schrodinger equation is merely possible and approxi- 
mative methods should inevitably be taken into account. 
The ground state wavcfunction of He-atom is node-less 
and the Pauli exclusion principle does not play a signifi- 
cant role. For this reason the above simple DMC method 
is reasonable. In table VI we have reported energy val- 
ues of He- atom in several ellipsoidal geometries. All used 
parameters such as a, b are defined similar to the H- 
atom case. By increasing a and b the ground state en- 
ergy asymptotically tends to the lowest state energy of 
free He-atom i.e. —2.903 au. In small radius the ener- 
gies have larger values, because electrons are closer to 
each other and to the nucleus. For higher values of Vq, 
two electrons are confined only inside the box and inter- 
act via Columbic forces, therefore the energy should be 
increased. The elimination of the electron-electron inter- 
action decreases the ground state energy in comparison 
to the case where this interaction is on. 

In the absence of nucleus, systems with the larger val- 



ues of a and b behave as a 2-dimensional confined electron 
gas which are vastly studied in quantum dots (li^] . 



VI. SUMMARY 

In this paper we have evaluated the ground state ener- 
gies and wavefunctions of confined Hydrogen and Helium 
atoms by the method of Diffusion Monte Carlo. In partic- 
ular, we have studied spherical and ellipsoidal confining 
compartments. Ground state values crucially depend on 
the size of confining compartments. In larger compart- 
ments, the results asymptotically approach to the values 
of free atoms. The effect of displacing the nucleus from 
the geometrical centre of the confining region has also 
been investigated. This effect is of theoretical as well as 
practical interest. We have shown that the ground state 
energies for both H and He atoms inside the ellipsoidal 
box significantly depend on the position of nucleus and 
the penetrability of boundaries. The merit of DMC is 
that it simply can treat an arbitrary shape of bound- 
aries. Our future work will be devoted to find a reason- 
able model for confined atoms within a carbon nano-tube 
and correspondingly applying QMC methods on such a 
systems with a cylindrical confinement. Work along this 
line is in progress and will be reported elsewhere. 
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a 


Vo Fi Y2 


Vo Fi Y2 Yz 


Vo Fi Y2 Fs 


Vo Fi Fa Fa 


Vo Fi Fa Fa 


Vo Fi Fa Fa 


1.0 
2.0 
3.0 
4.0 

5.0 


0.0 0.01 2.0 
0.0 0.04 2.23 
0.0 0.0 0.02 2.07 
0.0 0.42 2.61 
0.01 0.44 2.59 


0.0 0.22 3.19 
0.0 0.24 2.89 
1.0 0.03 0.0 2.14 
0.01 0.64 2.87 
0.01 0.74 2.90 


0.08 0.93 4.99 
0.01 0.47 3.41 
2.0 0.01 0.3 2.64 
0.01 0.31 2.54 
0.01 0.5 2.66 


0.01 0.9 4.68 
0.01 0.52 3.35 
3.0 0.01 0.24 2.57 
0.01 0.16 2.23 
0.01 0.48 2.51 


0.06 0.14 3.73 
0.01 0.62 3.31 
4.0 0.09 0.51 2.89 
0.03 0.61 2.75 
0.01 0.11 2.24 


0.02 0.95 4.93 
0.08 0.54 3.39 
5.0 0.01 0.83 3.31 
0.02 0.23 2.50 
0.01 0.41 2.61 



TABLE I: Values of the fitted parameters needed for trial wavefunctions in spherical penetrable box. 



a=l 


a=2 


a=3 


a=4 


a=5 


d A B 


d A B 


d A B 


d A B 


d A B 


0.1 0.035 0.498 
0.2 0.072 0.490 
0.3 0.106 0.478 
0.4 0.134 0.464 
0.5 0.152 0.447 


0.2 0.009 0.249 
0.4 0.018 0.245 
0.6 0.026 0.239 
0.8 0.033 0.232 
1.0 0.038 0.224 


0.3 0.004 0.166 
0.6 0.008 0.163 
0.9 0.012 0.160 
1.2 0.015 0.155 
1.5 0.017 0.149 


0.4 0.002 0.124 
0.8 0.005 0.123 
1.2 0.007 0.119 
1.6 0.008 0.116 
2.0 0.010 0.112 


0.5 0.001 0.099 
1.0 0.003 0.098 
1.5 0.004 0.096 
2.0 0.005 0.093 
2.5 0.006 0.089 



TABLE IL Values of coefficients A and B in sec. IV.B 




FIG. 1: (Color online) The electronic ground state wavcfunction of H-atom confined in a sphere for two values of Vq. a = 1 and 4.0 
(a) and (b) respectively. In both figures d = 0. 
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FIG. 2: (Color online) (a) Ground state energy of H-atom versus penetrability in an spherical compartment for various radii when d = 0. 
(b) Ground state energy of confined H-atom for different penetrability coefficients versus off-centre distances d when a = 1.0. (c) Ground 
state energy of a H-atom confined in an impenetrable spherical compartment for various values of a versus off-centre distance d for VO 
=00. 




FIG. 3: (Color online) Electronic distribution of H-atom inside a spherical compartment for Vq = and fixed a = 1 when d = 0. 




FIG. 4: (Color online) Electronic distribution of H-atom inside a spherical compartment for Vq = 5 and fixed a = 1 when d 
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FIG. 5: (Color online) Electronic distribution of H-atom inside a spherical compartment for Vq = /infty and fixed a = 



1 wiien d = 0. 
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a 


d 


Vo = 0.0 


Vo = 5.0 


Vo = oo 


1.00 


0.00 


-0.1250(7) 


0.822(7) 


2.1760(6) 




0.10 


-0.122(8) 


0.878(8) 


2.5070(7) 




0.20 


-0.116(9) 


0.918(4) 


2.5970(8) 




0.30 


-0.102(5) 


0.974(6) 


2.6210(6) 




0.40 


-0.082(9) 


1.040(6) 


2.711(5) 




0.50 


-0.059(7) 


1.1450(6) 


2.863(1) 


2.00 


0.00 


-0.430(3) 


-0.182(1) 


-0.1250(6) 




0.20 


-0.4280(7) 


-0.163(3) 


-0.027(8) 




0.40 


-0.414(2) 


-0.1500(8) 


-0.043(7) 




0.60 


-0.399(5) 


-0.113(6) 


-0.0260(3) 




0.80 


-0.394(2) 


-0.028(1) 


0.021(1) 




1.00 


-0.383(2) 


0.0600(8) 


0.110(2) 


3.00 


0.00 


-0.489(6) 


-0.4510(6) 


-0.425(7) 




0.30 


-0.4840(6) 


-0.446(4) 


-0.399(8) 




0.60 


-0.479(8) 


-0.4370(5) 


-0.385(6) 




0.90 


-0.467(4) 


-0.428(4) 


-0.384(3) 




1.20 


-0.463(8) 


-0.404(4) 


-0.371(3) 




1.50 


-0.4551(4) 


-0.351(9) 


-0.296(5) 


4.00 


0.00 


-0.496(3) 


-0.489(3) 


-0.481(9) 




0.40 


-0.493(7) 


-0.475(9) 


-0.470(3) 




0.80 


-0.491(8) 


-0.449(1) 


-0.465(2) 




1.20 


-0.489(9) 


-0.470(7) 


-0.459(4) 




1.60 


-0.480(7) 


-0.408(1) 


-0.463(5) 




2.00 


-0.470(4) 


-0.391(4) 


-0.435(4) 


5.00 


0.00 


-0.499(9) 


-0.498(6) 


-0.492(2) 




0.50 


-0.494(7) 


-0.490(2) 


-0.488(1) 




1.00 


-0.485(2) 


-0.479(6) 


-0.494(7) 




1.50 


-0.491(8) 


-0.478(8) 


-0.486(4) 




2.00 


-0.490(6) 


-0.469(4) 


-0.490(8) 




2.50 


-0.499(7) 


-0.471(5) 


-0.471(6) 



TABLE III: The ground state energy of confined H-. 
d. 



■atom for various values of penetrability coefficients and off-center distances 
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a 


d 


Vo = 0.0 


Vo = 5.0 


Vo = CO 


1.00 


0.00 


-2.2290(7) 


-0.910(5) 


1.012(4) 




0.10 


-2.228(8) 


-0.949(6) 


1.105(6) 




0.20 


-2.269(2) 


-0.965(2) 


1.175(9) 




0.30 


-2.241(8) 


-0.921(4) 


1.477(5) 




0.40 


-2.250(8) 


-0.941(6) 


1.926(2) 




0.50 


-2.247(4) 


-0.916(8) 


2.602(7) 


2.00 


0.00 


-2.757(5) 


-2.652(6) 


-2.603(6) 




0.20 


-2.774(7) 


-2.619(9) 


-2.496(3) 




0.40 


-2.747(7) 


-2.615(2) 


-2.448(7) 




0.60 


-2.723(2) 


-2.6300(8) 


-2.382(6) 




0.80 


-2.759(7) 


-2.642(8) 


-2.247(1) 




1.00 


-2.781(3) 


-2.622(9) 


-2.008(8) 


3.00 


0.00 


-2.831(4) 


-2.776(1) 


-2.871(8) 




0.30 


-2.797(8) 


-2.791(6) 


-2.754(6) 




0.60 


-2.801(1) 


-2.805(7) 


-2.749(4) 




0.90 


-2.781(9) 


-2.7820(3) 


-2.735(4) 




1.20 


-2.780(2) 


-2.757(2) 


-2.703(1) 




1.50 


-2.809(7) 


-2.806(6) 


-2.607(2) 


4.00 


0.00 


-2.816(5) 


-2.803(1) 


-2.892(6) 




0.40 


-2.771(6) 


-2.804(9) 


-2.800(1) 




0.80 


-2.806(8) 


-2.830(1) 


-2.795(8) 




1.20 


-2.7980(8) 


-2.787(4) 


-2.781(8) 




1.60 


-2.814(1) 


-2.792(9) 


-2.778(1) 




2.00 


-2.801(4) 


-2.789(9) 


-2.7380(4) 


5.00 


0.00 


-2.812(8) 


-2.798(2) 


-2.902(7) 




0.50 


-2.815(1) 


-2.806(1) 


-2.8710(8) 




1.00 


-2.8280(7) 


-2.776(2) 


-2.820(9) 




1.50 


-2.765(4) 


-2.7750(5) 


-2.797(5) 




2.00 


-2.797(9) 


-2.793(7) 


-2.812(8) 




2.50 


-2.821(6) 


-2.791(9) 


-2.800(4) 



TABLE IV: The ground state energy of confined i?e-atom for various values of penetrability coefficients and off-centre distances 
d. 



a 


b 


Vo = 0.0 


Vo = 5.0 


Vo — oo 


b 


Vo = 0.0 


Vo = 5.0 


Vo = oo 


1.00 


1.10 


-0.1320(2) 


0.819(6) 


1.832(1) 


1.50 


-0.179(2) 


0.605(9) 


1.460(7) 


1.50 


1.65 


-0.352(1) 


-0.013(1) 


0.269(8) 


2.25 


-0.367(9) 


-0.095(9) 


0.133(8) 


2.00 


2.20 


-0.428(6) 


-0.288(8) 


-0.184(3) 


3.00 


-0.443(7) 


-0.338(4) 


-0.257(8) 


2.50 


2.75 


-0.466(7) 


-0.415(1) 


-0.374(1) 


3.75 


-0.476(9) 


-0.428(7) 


-0.397(2) 


3.00 


3.30 


-0.482(1) 


-0.450(5) 


-0.434(3) 


4.50 


-0.485(4) 


-0.466(1) 


-0.452(2) 


3.50 


3.85 


-0.490(6) 


-0.473(4) 


-0.477(4) 


5.25 


-0.489(3) 


-0.480(2) 


-0.481(2) 


4.00 


4.40 


-0.490(6) 


-0.484(5) 


-0.492(3) 


6.00 


-0.492(3) 


-0.492(3) 


-0.494(6) 


4.50 


4.95 


-0.491(3) 


-0.490(6) 


-0.489(2) 


6.75 


-0.497(2) 


-0.493(1) 


-0.4940(6) 


5.00 


5.50 


-0.490(9) 


-0.490(2) 


-0.493(7) 


7.50 


-0.486(5) 


-0.497(8) 


-0.500(8) 



TABLE V: The ground state energy of H-atom confined in an ellipsoidal box, with b = 1.1a and b = 1.5a. 



a 


b 


Vo = 0.0 


Vo = 5.0 


Vo = oo 


b 


Vo = 0.0 


Vo = 5.0 


V() = oo 


1.00 


1.10 


-1.6890(3) 


-1.101(7) 


0.203(4) 


1.50 


-1.923(3) 


-1.429(3) 


-0.371(2) 


1.50 


1.65 


-2.544(6) 


-2.382(3) 


-2.115(6) 


2.25 


-2.611(9) 


-2.502(7) 


-2.274(3) 


2.00 


2.20 


-2.795(2) 


-2.750(4) 


-2.644(3) 


3.00 


-2.780(6) 


-2.772(8) 


-2.712(8) 


2.50 


2.75 


-2.881(4) 


-2.8350(9) 


-2.813(3) 


3.75 


-2.851(9) 


-2.851(6) 


-2.820(7) 


3.00 


3.30 


-2.859(7) 


-2.874(5) 


-2.880(3) 


4.50 


-2.904(3) 


-2.869(3) 


-2.8780(6) 


3.50 


3.85 


-2.868(3) 


-2.888(7) 


-2.908(2) 


5.25 


-2.885(2) 


-2.924(8) 


-2.866(8) 


4.00 


4.40 


-2.877(8) 


-2.894(4) 


-2.865(7) 


6.00 


-2.8900(6) 


-2.881(8) 


-2.885(9) 


4.50 


4.95 


-2.8480(1) 


-2.895(3) 


-2.889(2) 


6.75 


-2.883(9) 


-2.864(9) 


-2.860(5) 


5.00 


5.50 


-2.887(2) 


-2.9150(8) 


-2.899(2) 


7.50 


-2.889(4) 


-2.870(3) 


-2.869(3) 



TABLE VI: The ground state energy of He-atom confined in an ellipsoidal box, with b = 1.1a and b = 1.5a. 



